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We analyze spin transport in insulating antiferromagnets described by the XXZ Heisenberg model 
in two and three dimensions. Spin currents can be generated by a magnetic-field gradient or, 
in systems with spin-orbit coupling, perpendicular to a time-dependent electric field. The Kubo 
formula for the longitudinal spin conductivity is derived analogously to the Kubo formula for the 
optical conductivity of electronic systems. The spin conductivity is calculated within interacting 
spin-wave theory. In the Ising regime, the XXZ magnet is a spin insulator. For the isotropic 
Heisenberg model, the dimensionality of the system plays a crucial role: In d = 3 the regular part of 
the spin conductivity vanishes linearly in the zero frequency limit, whereas in d = 2 it approaches 
a finite zero frequency value. 
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I. INTRODUCTION 

The challenge of spintronics research is to exploit the 
spin degree of freedom as an additional tool in novel elec- 
tronic devices 0, 0, H, 0, [H . This task demands to ex- 
plore the basic physical principles underlying the genera- 
tion and decay of spin-polarized charge currents. In this 
context, metallic and semiconducting devices have been 
examined both theoretically and experimentally. Simul- 
taneously, the synthesis of quasi-one-dimensional corre- 
lated insulators such as Sr 2 Cu03 Q and CuGeC>3 has 
renewed the general interest in the spin and thermal 
transport properties of low-dimensional quantum spin 
systems. Also fundamental questions have been raised, 
how the integrability of a system influences its spin 
g i, 01, 03, El and thermal 0, 0, Q 01 conductiv- 
ities. In two-dimensional high-mobility electron systems 
with Rashba spin-orbit coupling charge currents are nec- 
essarily accompanied by spin currents [l6l |. These spin 
currents flow perpendicular to the charge current direc- 
tion and therefore lead to an intrinsic spin Hall effect 

MM 

The spin conductivity of Heisenberg chains has previ- 
ously been computed within linear-response theory by 
adopting an analogy to the Kubo formula for charge 
transport [l9j . The low- frequency behavior of the optical 
conductivity <j{lo) = + ia"(uj) provides a transpar- 

ent scheme to distinguish the charge transport properties 
of ideal conductors, insulators, and non-ideal conductors 
[l9L [20| . Decomposing the real part of the longitudinal 
conductivity as cr'((d) = D8(oj) + a ros (ui), a finite Drudc 
weight or charge stiffness |2l| D > is the character- 
istic of ideal conductors. A similar classification scheme 
can be carried over to spin transport and the spin con- 
ductivity in order to distinguish between spin conductors, 
spin insulators or even spin superfluids [23 . |23| . 

The optical conductivity is conveniently derived as the 
current response to a time-dependent electromagnetic 
vector potential. Similarly, for the spin current response 
the concept of a fictitious spin vector potential can be 



introduced, which is related to a twist in the direction 
of the spin quantization axis [24j. However, the phys- 



ical realization of this perturbation and its relation to 
externally applied magnetic or electric fields is not obvi- 
ous. Indeed, spin currents flow in response to a magnetic- 
field gradient. In this case, the analogy to the generation 
of electric currents by a potential gradient is straight- 
forwardly established for one-dimensional systems by a 
Jordan- Wigncr transformation [25|, [2(| [27} ■ This trans- 
formation maps spin-i operators to creation and anni- 
hilation operators of spinlcss fermions; a magnetic-field 
gradient for the spins thereby translates into a potential 
gradient for the fermions. 

In this paper, we do not appeal to the analogy to the 
charge current response. Rather, we derive the Kubo for- 
mula for the spin conductivity of XXZ Heisenberg mag- 
nets directly (Sec. |n|. In particular, a magnetic-field 
gradient or, in the presence of spin-orbit coupling [28j |. 
a time-dependent electric field can be used to drive a 
spin current (Sec. IHIj) . The longitudinal spin conductiv- 
ity of the Heisenberg antiferromagnet in two and three 
dimensions is computed using spin- wave theory (Sec. IIVP 
for both the non-interacting-magnon approximation (Sec. 
IV B[) and the ladder approximation for repeated two- 
magnon scattering processes (Sec. IV Cjl . In particular, 
the low-frequency behavior of the spin conductivity will 
be analyzed in detail revealing distinct differences be- 
tween the spin transport properties of 2D and 3D anti- 
ferromagnets. 



II. KUBO FORMULA 

Specifically, consider the antiferromagnetic XXZ 
Heisenberg model (HAFM) 



(1) 



with nearest-neighbor exchange coupling J > 0, an 
anisotropy parameter A, and local spin operators Si of 



2 



length S. The sum over extends over the nearest- 
neighbor bonds on a d-dimensional hypercubic lattice 
with lattice constant a = 1 and N sites. An external 
space- and time-dependent magnetic field B z (l,t) cou- 
ples to the spin system via the Zeeman energy. The 
time-dependent Hamiltonian thus reads 



H(t)=H-^S z (l)h z (l,t) : 



(2) 




FIG. 1: (Color online) Setup for a spin current generated by 
a magnetic-field B z gradient along the x direction. 



where h(l,t) = h z (l,t)e z = g^sB z (l,t)e z . 

The spin current density operator ji^j for the magne- 
tization transport from site i to site j is defined via the 
continuity equation 



dtSf 



(3) 



Here, J2j h^j 1S the lattice divergence of the local spin 
current density at site i. The operator j x {l) = ji^i+ x 
thus follows from Heisenberg's equation of motion S z = 
i[H, Sf] and Eq. © as 



Jx(0 



Ji 



l+x 



(4) 



where I + x is the nearest-neighbor site of site I in the 
positive x direction. 

In our subsequent calculation we assume long-range 
antifcrromagnctic order oriented along the z direction 
and thus restrict the analysis to a scalar spin current op- 
erator for the magnetization transport. The definition of 
a proper spin current vector operator Jj_ >j is rather sub- 
tle, as discussed in detail in Ref. [2j|. In fact, if (Si) and 
(Sj) are not collinear, only the projection of Ii—>j onto 
the plane spanned by the local order parameters (Si) and 
(Sj) may be interpreted as a physical transport current. 
If (Si) and (Sj) are collinear, however, the magnetization 
transport is indeed correctly described by the scalar cur- 
rent density operator defined in Eq. ((4]). We note that 
the proper definition of a spin current operator is also 
a controversial issue in the context of spin transport in 
semiconductors with spin-orbit coupling [30], but these 
issues are of no concern for the purpose of our analysis. 

The linear spin current response to an external 
magnetic-field gradient is 



(jx{q,u)) = x 3 s{q,u)h z (q,u), 



(5) 



where the Fourier transforms in time and space are de- 
fined as 

/oo 
dte ^t- q -l) A{ht) (6) 

In Eq. ([5]) we have introduced the dynamic susceptibility 

XjS (q,u) = — j a dte^ +io+ »(\j x (<l,t),S z (-q,0)})- 

(7) 



Using the spatial Fourier transform of the continuity 
equation 



S z (q,t)+iq-j(q,t) = 0, 



(8) 



and assuming a magnetic-field gradient of B z only along 
the x direction, \jS is transformed by partial integration 
as 



1 



HN i(u> + i0+) 



([Mq),s z (-q)}) 



(9) 



+ / &te^ +i0+ »(\j x (q,t)Aq x j x (-q,Q)]) 
Jo 



The response formula therefore becomes 



(jx(q,u)) 



-K r . 



K x (q,w). 



i(cj + i0+) 



iq x h z (q,Lo), (10) 



where we have introduced the spin-flip part of the ex- 
change interaction along the x direction 



(K x ) 



^£j<S+S- x + S r S+ x >, (11) 



2hN 



and the longitudinal retarded current-current correlation 
function, 

A xx (q,u;) = — I dte^ +io+ »(\j x (q,t),j x (-q,0)]). 

(12) 

For a closer analogy to charge transport, we consider 
the transport of magnetization m; = g^sSf instead of 
(dimcnsionlcss) spin and define the magnetization cur- 
rent operator j m x = g^Bjx- The longitudinal spin con- 
ductivity a xx (uj) is then defined as the linear magneti- 
zation current response to a long-wavelength (q —> 0), 
frequency dependent magnetic-field gradient (Fig. [1|. 
Putting h z = gnsB z , the relation 



jm, x (q,v) = (T xx (q,u))iq x B z (q,ui) 



(13) 



thus yields the Kubo formula for the spin conductivity 
in the long- wavelength limit, 



2 \~ K X ) - A xx (q = 0,u) 
i(oj + i0+) 



(14) 



This result indeed establishes the perfect analogy to the 
formula for the optical conductivity of interacting lattice 
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electrons [l9[. The real part of the spin conductivity is 
given by 



(15) 



with the spin stiffness [22|, |23[ or "spin Drude weight" 
D S , 



D f 



{ gi i B f{(-K x ) - A' xx (q = 0,w - 0)), (16) 



and the regular part 



(17) 



The spin conductivity also fulfills the "f-sum rule" [3l| 



duj °'xx{u) = {-K x ), 



(18) 



which can be derived by using the Kubo formula and the 
Kramers-Kronig relations for A xx . Note that the integral 
in Eq. (fT5|) contains one half of the possible spin Drude 
weight peak at zero frequency. 

The structure of the spin conductivity formula emerges 
from the straightforward calculation for the linear cur- 
rent response to an external magnetic-field gradient with- 
out introducing Pcierls-like phase factors with a fictitious 
spin vector potential 0, [32j|. However, in an external 
electric field a moving magnetic dipole does acquire an 
Aharonov-Casher phase [33f . Katsura et al. [28| pointed 
out that Aharonov-Casher phase factors and a corre- 
sponding "vector potential" can be introduced on the 
basis of the Dzyaloshinskii-Moriya interaction [34|, [3^ . 
which may be induced in an external electric field due to 
spin-orbit coupling. As we show in the next section, this 
alternative approach leads to the same result for the spin 
conductivity, Eq. (CO 



III. RESPONSE TO AN ELECTRIC FIELD 

Spin currents are also generated by applying a time- 
dependent electric field in the presence of spin-orbit cou- 
pling. In low-symmetry crystals with localized spin mo- 
ments spin-orbit coupling, parametrized by a coupling 
constant A so , leads to the Dzyaloshinskii-Moriya (DzM) 
antisymmetric exchange interaction [H, [H| 



zM 



(id) 



(Si x S 3 ). 



(19) 



In high-symmetry crystals an inversion symmetry break- 
ing external electric field E induces a DzM vector cx 



E x eij [36 



where e^j denotes the unit vector connect- 



ing the neighboring sites i and j. Specifically, for a field 
E = (0, E y ,0), the DzM vector on the bonds in the x-y 
plane takes the form 



D 



aE x e a 



-aE y e z , 



D 



i,i-\-y 



o. 



(20) 



where a cx A so . With (£< x SjY = ±{S+SJ~ - SrS+), 
the DzM Hamiltonian can therefore be rewritten as 



(21) 



where j x (l) is the spin current density operator defined 
m Eq. gl) and 



A x (l,t) 



2aE y (l,t) 
J ' 



(22) 



The electric field therefore acts as a "spin vector poten- 
tial" in formal analogy to the electromagnetic vector po- 
tential in charge transport. 

If the DzM exchange interaction is added to the XXZ 
spin Hamiltonian, the Heisenberg equation of motion 
yields an additional contribution to the spin current op- 
erator, 



•DzM 

Ji— >j 



1 



(D Jl -D lJ )-(S+Sr+SrS+), (23) 
and the total magnetization current operator is given by 
jm,xQ; A) = gfi B jx(l;A = 0) + {gfx B ) 2 K x (l)A x {l). (24) 

As for the charge current response to an electromagnetic 
vector potential [l9| , the linear spin current response to 
the spin vector potential in the long-wavelength limit be- 
comes 

(j m ,x(°,u)) = (9Vb) 2 ((K x ) +A xx (0,lj))A x (0,u). 

. (25) 

The spin conductivity is therefore obtained by identifying 
the time derivative of the spin vector potential, i(u> + 
i0 + )A x (O, uj), as the driving force for the spin current. 
As a consequence, in the presence of spin-orbit coupling, 
a spin current can also be driven perpendicular to a time- 
dependent electric field (see Fig.[2]) . Physically, the origin 
of this phenomenon is contained in the Maxwell equation 



„ 4tt . IdE 
rot B = — j c H — , 

c c ot 



(26) 



in the absence of charge currents, j c = 0. Assuming a 
magnetic field in z direction, the y component of Eq. (|26p 
reduces to 



d x B z = —d t E v = —^-d t A x . 
c v 2ac 



(27) 



Eq. (|27|) identifies the time-dependent electric field and 
the magnetic field gradient as the same driving forces for 
the spin current. 



IV. SPIN- WAVE THEORY FOR THE 
ANTIFERROMAGNETIC XXZ MODEL 

In our subsequent analysis of the dynamic spin cur- 
rent correlation function A xx (0, to) we employ interacting 
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■d/dt E v 




FIG. 2: (Color online) Setup for a transverse spin current 
driven by a time-dependent electric field in the presence of 
spin-orbit coupling. 



For a d-dimcnsional hypcrcubic lattice the spin-wave dis- 
persion is 



tM k = 2dAJSa(S)e k , e k = y/l - 7 £/A 2 (35) 

with 7^ = 2 Sq=i cos(fc Q ). The magnon- vacuum energy 
is E = —NAJS 2 a 2 (S)d. The Oguchi correction factor 



(36) 



spin-wave theory using the Dyson-Maleev (DM) transfor- 
mation [H, [H, [4(| for the antiferromagnetic Heisenberg 
model on bipartite lattices. In the DM representation 
the spin operators are replaced by bosonic operators ac- 
cording to 



— S (1^ (Xi , 



(28) 



for the up-spin sublattice A and by 

= -S + bfa, (29) 

for the down-spin sublattice B. Due to the bosonic com- 
mutation relations for the a and b operators, the spin 
algebra is preserved. 

The spatial Fourier transformation of the bosonic op- 
erators is conveniently written as (4l| 



where the momentum k is restricted to the magnetic Bril- 
louin zone (MBZ). For the 2D square lattice MBZ = {k : 
I Mb | + \ky\ < tt} ■ Inserting the DM representation into 
the XXZ Heisenberg model, the quadratic part of the 
resulting bosonic Hamiltonian is diagonalizcd by the Bo- 
goliubov transformation 

a k = u k a k + v k Pl, b k = u k (3 k + v k a\. (31) 
The coefficients u k and v k are given by 



l + £k . . /1-Sfe 

u k = \l—^i — , Wfc = -sgn(7feJ 



(32) 



2e k V 2gfc 

The DM transformed Hamiltonian reads 

Hdm = E + H + Vdm (33) 
with the diagonalizcd quadratic part 

Ho = ^fi^(4ajc + ^A). (34) 

k 



arises from the normal-ordering of quartic terms [42J in 
the interaction part Vdm- In d = 2, 



1 



3F2 



111 1 



(37) 



with the generalized hypergeometric function 3F2 [43[ | . 
Specifically, for d = 2, A = 1, and S = |, the Oguchi 
correction factor is given by 



a(i) 



r(i) 4 



r(|) 4 



8- 



1.157947, 



(38) 



where T denotes the gamma function. Finally, the 
normal-ordered quartic interaction part of the Dyson- 
Maleev Hamiltonian (l33l) reads 



(39) 



Vdm = -A J- ]T <5 G (1 + 2-3-4) 



(1234) 



V [1) a\a\a z a 4 + V (2) a\l3 2 a 3 ai + V A(3) a t 1 a^/3|a 4 
+y (4) a t 1 a 3 /5l/?2 + V^(3\a^x + V^Malfii 



The interaction vertices i = 1, . . . ,9 depend on the 
wave vectors k\ , . . . , k±, abbreviated by 1, . . . , 4, and are 
explicitly given in Ref. [33]. The Kronccker delta, <5g(1 + 
2 — 3 — 4), ensures momentum conservation to within a 
reciprocal lattice vector of the MBZ. 

We briefly comment on two difficulties with the Dyson- 
Maleev transformation. Firstly, two constraints, a|a» < 
2S and bjbj < 2S, arc required for the Fock space of 
bosonic excitations in order to avoid unphysical spin ex- 
citations. The spin- wave analysis without constraints 
can nevertheless be quantitatively justified noting that 
(ajai) sa 0.197 in d = 2 at T = for the S = \ isotropic 
Heisenberg model, calculated in linear spin-wave theory 
(LSW), where Vdm is neglected. This result implies for 
the sublattice magnetization (jtia) = \ — (fljoii) ~ 0.303 
(A = 1), supporting a posteriori the validity of spin- wave 
theory even for S = ^. 

Secondly, the Dyson-Maleev transformation is not Her- 
mitian, since S + 7^ (S~Y when the spin operators are 
expressed in terms of bosonic operators. As a conse- 
quence, the quartic part of the transformed Hamiltonian 



5 



is non-Hermitian, too. However, its non-Hermiticity does 
not affect our calculations since we will take into account 
the only quantitatively relevant Hcrmitian V^ 4 ^-term (see 
Sec. ES| . 



V. EVALUATION OF THE SPIN 
CONDUCTIVITY 

A. Basic propagators and correlation functions 

We proceed by defining the time-ordered magnon prop- 
agators (see Refs. [13, EJ for more details), 



G aa (fc, t) 



-i(0|Ta fc (*)4(0)|0), 



(40) 



G^(k,t) = -i<0|T/3£(i)/3 fe (0)|0>. (41) 

The bare Fourier-transformed propagators in the absence 
of interactions are 



G^(k,uj) = 



1 



w - fi fc + i0+ 



-1 



w + fife - i0+ 
. . (42) 

The regular part of the longitudinal spin conductivity 



(,9MB ) 2 

UJ UJ 



2 A^(q = 0,w) 



G» (43) 



is determined by the imaginary part of the Fourier trans- 
formed time-ordered spin current correlation function 



G 3 (t) 



hN 



(0\T Jx (t) Jx (0)\0) 



(44) 



in the ground state |0) with the spin current operator 
j x = J2i jx(l) (see Eq. ([3])). j x is transformed by the same 
steps which are used for the Dyson-Maleev transforma- 
tion of the HAFM Hamiltonian. The result is j x = jo+ji, 
where 



Jo 



and 



2JSa(S)J2 s Hkx) 



(45) 



7fc 



.71 



J ^ Z fc(l + 2-3-4)sin(/ci, x ) 



(46) 



(1234) 



+j (4) a t 1 a 3 /3l/3 2 + J [b) P\a^ 2 Pi + M<£Pi 



where, for brevity, we have omitted here the explicit ex- 
pressions of the spin current vertices j'W; for complete- 
ness, they are listed in the Appendix. The quadratic part 



jo is of order S 1 , but it also contains Oguchi corrections 
of order S° arising from the normal-ordering of quartic 
terms in j\ . 

The calculation of Gj (t) involves the following four cor- 
relation functions: 

(i) A two-magnon correlation function 



Gf(t) = -i<0|Tjo(t)io(0)|0), 

(ii) a four-magnon correlation function 

G?(t) = -i{0\Tj 1 (t)j 1 (0)\0), 

(iii) and two cross-correlation functions 

Gf(t) = -i(0|Tjo(t)ji(0)|0), 
Gf (t) = -i<0|r7i(*)jo(0)|0>. 



(47) 

(48) 

(49) 
(50) 



The regular part of the spin conductivity (|43|) is then 
obtained with 

G» = \ [G"° (uj) + Gf{u)) + Gf(u) + . 

(51) 

Before we continue with the calculation of Gj(uj), we 
discuss the selection rules for the matrix elements of the 
spin current operator, which provide insight into the rel- 
evant physical processes for spin transport in antifcrro- 
magnets. The Lehmann representation 



1 E n ^E 

G » =iiE l(0bx|n}| 2 



(52) 



uj - (E n - E ) + i0+ uj + (E n - E )- i0+ 

shows that the spin current correlation function is deter- 
mined by the matrix elements (0|ja;|n), where |n) denotes 
the excited states with energy E n . The ground state 
and the relevant excited states must therefore have (i) 
the same total momentum and (ii) the same spin quan- 
tum numbers, but (iii) opposite parity in order to have 
nonvanishing matrix elements (0|ja;|n). These conditions 
imply that only multi-magnon excitations with vanish- 
ing total momentum contribute to the spin conductivity; 
one-magnon excitations are forbidden by both (i) and 
(ii). Condition (iii) is reflected in the sm{k x ) vertex in 
Eq. (|46|) . which selects the x direction for the spin cur- 
rent. Below we will focus on the two-magnon contribu- 
tion to the spin conductivity. Apart from the sm(k x ) 
vertex function, our calculation proceeds analogously to 
the two-magnon analysis of Raman scattering in anti- 
ferromagnets [37j]. In fact, the selection rules for spin 
transport are the same as for the two-magnon Raman 
scattering intensity in B\ g scattering geometry 

In the analysis of the spin current correlation function, 
Eq. (|5"T]) , we focus on G°° because it provides the leading 
(S 2 ) order contribution and also the dominant S 1 correc- 
tions for the spin conductivity. The product jo(t)jo(0) 
contains 16 terms. However, following the arguments of 
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Canali and Girvin in Ref. [37f , the contributions of 14 of 
these terms are negligibly small due to two arguments: 
Firstly, the 12 terms containing prefactors ^ k /Ae k are 
small because the absolute value of this factor is small 
near the MBZ boundary, where the free magnon density 
of states has a van- Hove singularity. Specifically, the four 
regions of the MBZ boundary where | sin(fec)| ~ 1 provide 
the quantitatively relevant contributions. Secondly, only 
two of the four remaining terms are nonzero if magnon 
interactions are neglected. 

The two remaining dominant terms lead to 



where 



(2JSa{S)y 
HN 



^sin(fc : )sin(fc,;) n 



k.k' 



£k£k> 



(53) 



(54) 



The two terms which are involved in Eq. [53] are thus 
calculated from the same two-magnon Green function, 

Hfcfc'C*) = -i(0|Ta fe (i)/3 fe (t)4'(0)/3L(0)|0). (55) 

Its Fourier transform is expressed in terms of magnon 
propagators and a 3-point vertex function pij 

/oo 
——G aa (k, uj+uj')Gpf3{k, Lo')T kk > (lo, lo'). 
-oo 2lT 

(56) 

The vertex function (lo, lo') satisfies the Bethe- 
Salpeter equation 



kk 



3fcfc' 



i AJd ^ 
h N 4-' 



dwi 
~2^ 



xV fcf 1 fc 1 fc( W '- W l) G "a(fcl,^+^l) 

*Gpp(ki, uJi)T klk >(uj, uji). (57) 

Here, the 4-point vertex V%£ fci k (u)',LOi) is the sum of all 
the irreducible interaction parts [44[ ■ 

In terms of magnon propagators and the vertex func- 
tion, Gj(u) is therefore given by 



i{2JSa(S)Y 



sin(k x ) 



hN j_ OQ -Zn <-£ e k 
xG aa (k, lo + w')Gpp(k, w')T h (w, w'),(58) 
where the reduced vertex function T k is defined by 



T k (u, lo') 



J2 s -^r kk ,(^<) 



£fe< 



(59) 



Eqs. ([57]) and ([58]) define the integral equations which 
have to be solved in the calculation of the correlation 
function G~^(lo). We will proceed in two steps: In a 
first approximation, the interactions between magnons 
are omitted (Sec. IV Bp . Subsequently, interactions 
are treated within a ladder approximation for repeated 
magnon- magnon scattering processes fScc. IV C]) . 
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FIG. 3: (Color online) Real part of the longitudinal spin con- 
ductivity for non-interacting magnons in two dimensions for 
different values of the anisotropy parameter A. The dashed 
lines show the expansion around u> gap , Eq. (I64|l . 



B. Non-interacting magnons 

We start in a first step by neglecting magnon-magnon 
interactions, which amounts to the replacements G — > 
G (0) and T fc -► sm(k x )/e k , i.e. V aP = 0. The required 
complex contour integral in Eq. ([57]) then straightfor- 
wardly leads to the result 



+ _ (JSa(S)f \ - sin 2 (fc x ) 
G * [U}) - hN ^~^~\ 



- 20* 



(60) 



For convenience we introduce the dimcnsionlcss fre- 
quency lo = w/O ma x, where f2 max = 2dAJSa(S)/h is 
the maximum one-magnon energy, and for m £ {0,1,2} 
we define the functions 



2 sm 2 (k x ) 



1 



N 



uj-2e k + i0+ 



(61) 



The regular part (tu > 0) of the spin conductivity within 
LSW is then given by 



h (dA) 2 LO 



lm.£ {2) (L 



(62) 



(g^B ) 2 



(dA) 2 u N 



E 



sin (k x ) 

2 o(w - 2£ fc ). 



Eq. ([62]) directly reflects the spin current selection rules 
explained above: 8(uj — 2e k ) accounts for two magnon 
excitations at energy e k . Momentum conservation and 
spin conservation arc fulfilled by a combination of a and 
f3 magnons, which carry opposite spin (S z = +1 or —1) 
and momentum (fc and —k). In a fully polarized Heisen- 
bcrg ferromagnet, there is only one magnon species and 
hence the spin conductivity vanishes @. 
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FIG. 4: (Color online) Regular part of the longitudinal spin 
conductivity within the ladder approximation for magnon- 
magnon interactions for S = | in d — 2. The dashed lines 
represent the linear expansions around w ga p, Eq. (JTTJ . 



The spin Drude weight (Eq. l[T5|l ) Ds = for any 
A > 1 within LSW. Furthermore, the regular part of the 
spin conductivity diverges at the maximum two-magnon 
energy (uj = 2) in dimensions d — 2 and d = 3 due to the 
neglect of magnon interactions. The LSW result for the 
spin conductivity in d = 2 is shown in Fig. [3] A special 
feature of the regular part of a' xx (uj) is its finite zero fre- 
quency limit at the isotropic point A = 1. An expansion 
near the magnon energy gap w gap = 2-^1 — 1/A 2 yields 
the leading contribution as 



{gfi B ) 2 /h 



A d di- 1 (5 2 



2¥ 



7T, d = 2 
1, d = 3 



(63) 



for cj > Wgap, with corrections of order (ui — w gap ) 
Specifically in d = 2. the result of the expansion is 



(64) 



for uj > lu„ 



The finite zero frequency 



value of the spin conductivity for A 

as 



1 is thus obtained 



</ xx (w - 0) = 



LSW 



h 8' 



(65) 



The expansion (|63[) for the spin conductivity within LSW 
in <i = 2 is also shown in Fig. [3l This expansion reveals 
that in d = 3, as opposed to the two-dimensional case, the 
spin conductivity for A = 1 vanishes linearly in the zero 
frequency limit within LSW (see Fig. [5]). The unphysical 
divergence at the upper edge of the LSW spectrum will 
be cured by including magnon-magnon interactions as 
discussed in the next subsection. 



C. Ladder approximation: two-magnon scattering 

For the calculation of C?1~(cj) magnon-magnon inter- 
actions are taken into account to lowest order [37j by 
approximating the 4-point vertex V a @ by its first-order 
irreducible interaction part, 



lja/3 / / \ _ T/ (4) 

v kk 1 k 1 k\ L0 > w iJ - v kk 1 k 1 k^ 
which is explicitly given by 



kk\k\k 



7fe-fei 



1 



£fc£fei 



1 - 



7fc7fci 
2A 2 £ fc e fel 



(66) 



(67) 



i.e., we neglect all the contributions to V al3 where two 
or more of the bare interactions V 1 - 1 ' are involved. Then 
the magnon propagators can again be replaced by the 
bare expressions gL°2 and g! ] in Eq. (|57j) since all the 
first order diagrams for the magnon self-energy vanish at 

t = o [13. 

The algebraic solution of the coupled integral equa- 
tions (|57p and (|5"8")) is based on the decoupling of the 
sums over k and k\ by means of the identities 

X] sin ( fc zhfcfffc = 0, (68) 

k 

^2sm(k x )j k ^ kl g k = sm (^i^) y^ s j n 2 (fc :i ,) gfci (69) 
k d k 

which hold for any function g k which has the symmetry 
of the hypercubic lattice. 

We proceed by repeatedly using Eqs. (f6"8"| and (|6"§|) to 
decouple Eqs. (|57|) and (|58|) leading to a ladder approxi- 
mation of the Bethe-Salpeter equation for the vertex part 
T k (uj,u)'). This has been demonstrated in Ref. (3?| for 
the vertex function (cos^) — cos{k y ))/2 of the Raman 
Big mode instead of the spin current vertex sin(fca;). By 
virtue of the identities Eqs. (|68|) and (|69|) the analogous 
analytical steps can be performed, and we obtain for the 
spin conductivity 



M = 



x Im 



{dA) 2 u 

£(2) „ K (£(1)£(D 



(70) 



1 + k(£(°) + £( 2 )) - k 2 (£W£W - £(o)£(2)) ' 



where = £^(u~>) and n' 1 = 2dSa(S). 

Figs. [4] and [5] show the S = | spin conductivity within 
the ladder approximation for d = 2 and d — 3, respec- 
tively. Scattering between magnons removes the diver- 
gence of the LSW result at ui = 2£! max . The frequencies 
for the maximal spin conductivity increase with increas- 
ing anisotropy parameter A. 

In contrast to the spin conductivity of the 3D Heisen- 
berg antiferromagnet, which vanishes at the isotropic 
point A = 1 for uj — > 0, the most notable feature of 
the regular part of the spin conductivity of the isotropic 
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FIG. 5: (Color online) Regular part of the longitudinal spin 
conductivity within the ladder approximation for S = \ in 
d — 3. The dashed lines show the LSW results. 



2D Heisenberg antiferromagnet remains its finite value 
in the zero frequency limit. For lu > w gap we find the 
following leading-order linear expansion in d = 2, 

^( w )^U cr 4(A)(l-^), (71) 

which is also included in Fig. [H A(A) is a numerical 
prefactor with A(l) = 1, and the zero frequency limit of 
the spin conductivity for A = 1 is given by 



coupling, perpendicular to a time-dependent electric 
field. We presented an explicit derivation of the Kubo for- 
mula for the spin conductivity for the antifcrromagnctic 
XXZ Heisenberg model and showed that the magnetiza- 
tion transport arises from two-magnon processes, which 
provide the dominant contribution to the spin conduc- 
tivity. 

In close analogy to the calculation of the B\ g two- 
magnon Raman light scattering intensity, the spin con- 
ductivity was evaluated using interacting spin-wave the- 
ory. The dimensionality of the model is important for 
the low-frequency behavior of the spin conductivity, es- 
pecially upon approaching the isotropic point (A = 1) 
from the gapped Ising regime (A > 1). In d = 3, the 
spin conductivity vanishes for A = 1 in the dc limit. In 
d = 2, however, the dc spin conductivity remains finite 
at the isotropic point. 

Experimentally the spin conductivity can be deter- 
mined by measurements of magnetization currents. This 
issue was discussed by Meier and Loss in Ref. [6j|, where 
possible experimental setups were proposed. In these se- 
tups the magnetization current is detected via the elec- 
tric field which it generates. Given the estimates in Ref. 
for the expected voltages in moderate magnetic field 
gradients, measurements of the spin conductivity indeed 
appear experimentally feasible. 

Acknowledgment: This work was supported by 
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► 0) 



ladder 



Yn 



h 8' 



1.856851 for S = -. 



(72) 
(73) 



APPENDIX: SPIN CURRENT VERTICES 



The renormalization factor Y a can be expressed in terms 
of the gamma function, 



Y„ 



1 — ban 



±h 2 K 2 

4°0 hu 



{l-\(b Q +b 2 )K-\{bl-b b 2 )K 2 y 



(74) 



32tt 



2-1 

7T 



47T 3 



16tt 



The inclusion of magnon-magnon interactions thus renor- 
malizes a* from its value for non-interacting magnons. 
The fact that <j' xx {uj) is finite in the limit w — > 0, how- 
ever, turns out to be a robust feature of the 2D Heisen- 
berg antiferromagnet. 



VI. SUMMARY AND CONCLUSIONS 



The spin current vertices involved in Eq. (I4T[) are given 



= 

/( 3 ) = 
;(4) . 



?' (5) = 
= 

,-CD = 

f( 8 ) = 

;(9) - 



U1M2U3U4 + V1V2V3VA + + U1U2V3V4, 

2(lti«2«3W4 + V1U2U3V4), 
2(uiU 3 U4M 2 + U1W3U4U2 + 

V1U3V4V2 + VIV3U4V2), 

2(v i u 3 u 2 vi + U4W3W2U1), 

V4V3V2V1 + U4M3M2U1 + V4V3U2U1 + U4U3W2W1, 

U1V2V3V4 + V1U2U3U4, 

V1U2U3U4 + U1V2V3W4, 

V4V3U2V1 + U4U3V2U1. 



In Heisenberg magnets, spin currents flow along a 
magnetic-field gradient or, in the presence of spin-orbit 



The coefficients m = Uki and Vi — with i = 1, . . . , 4 
are defined by Eq. (132]) . 
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